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1. Strangeness as a messanger 

Strangeness production in heavy-ion collisions (HICs) enjoys a persistent interest. It is ex¬ 
plained by several reasons. First, strangeness is a tag on a hadron, saying that this particular hadron 
was most probably produced in the course of collision and does not stem from the original nuclei. 
Second, strange quarks and anti-quarks prefer to reside in different hadronic species. A strange 
quark prefers to be in baryons: so there is two kaon K states for several baryonic states A, £ ,E, ££ 
Oppositely, the anti-strange quark can only be in mesons, like K and its heavier resonances. This 
difference leads to distinct interactions of strange and anti-strange subsystems with a non-strange 
baryonic environment. The strange subsystem couples more strongly than the anti-strange one. 
As argued in Ref. [1, 2, 3] this could leads to strangeness/anti-strangeness separation in baryon- 
rich matter formed in HIC at energies from SIS to middle RHIC energies. Third, strangeness is 
conserved in strong interactions. Therefore the strange and anti-strange particles are created in 
elementary processes only pairwise. This renders high thresholds for strange particle production, 
making yields of strange particles in HIC sensitive to possible in-medium effects. 

In spite of the common enthusiasm about strangeness in HIC, such studies are complicated by 
several factors. First of all, cross sections of reactions with the strangeness production are poorly 
known, especially close to the reaction threshold. Experimentally accessible are only reactions 
with the pp, np and np{n) entrance channel, as those currently studied, e.g., by HADES collabo¬ 
ration at GSI [4] and by ANKA, HIRES, and TOF collaborations at COSY [5]. Without reliable 
experimental information about reactions with the nn initial state we have to rely on the isospin 
symmetry relations, making ad hoc assumptions about internal mechanisms of the reactions [6]. 
Second, interactions of strange particles with non-strange ones and among themselves cannot be 
fully constrained. There are limited data in the channels with the total strangeness S = ±1, like a 
hyperon-nucleon and anti-kaon-nucleon scattering. Finally, the properties of strange particles and 
their interactions in dense and hot nuclear matter are far from being firmly established. 

Against all odds mentioned above a satisfactory description of the kaon and A productions in 
HICs has been reached in various transport [7, 8] and statistical [2, 9] approaches. 

There remains however several open issues with the production of particles containing several 
strange quarks, like E and Q. baryons and 0 mesons, which we would like to discuss in the present 
talk. 

2. Strange particles in medium 

Medium effects prove to be important for description of particle production in HIC at SIS en¬ 
ergies [10, 1 1]. Particularly, strange particle yields are strongly influenced by them [7, 8, 2, 12, 13]. 
For hyperons and nucleons the in-medium modification of the energy spectrum of particle a is often 
effectively parameterized in terms of scalar S a and vector V a potentials E a (p) = sjm * 2 + p 2 + V a , 
where the scalar potential enters the spectrum through the effective mass m* = m a +S a . Description 
in terms of the S a and V a potentials is typical for relativistic mean-field (RMF) models, cf. [7, 14]. 
The nucleon potentials determined in [14] are .S' ; y ~ —190 MeVp^/po and Vn — +130 MeVp^/po, 
here ps is the baryon density and po = 0.16/I ni ' is the nuclear saturation density. The same poten¬ 
tials could be also used for A: Sa — Sn, Va — Vn. 
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Figure 1: World of kaon-nucleon (upper panel) and anti-kaon-nucleon (lower panel) scattering. Vertical 
dashes show thresholds and fat dots correspond to poles. 


One usually relates the hyperon potentials to the nucleon ones, Vh = HhVn and Sh = PhSn- 
The parameter a# one choose according to the number of non-strange quarks in the hyperon, G!a = 
otv = 2az =2/3. The parameter /3 h is, then, chosen such that the optical potential of a hyperon in 
nuclear medium at saturation Uh = Sjj(po) — V#(po) agrees with the empirical information from 
the analysis of hypernuclei: U\ = —27 MeV [15], Uz = +24 MeV [16], and Uz = —14 MeV [17]. 

In the sector of strange mesons the situation is more controversial. In Fig. 1 we depict maps 
of the worlds of KN and KN scattering. The KN interaction looks structureless without any res¬ 
onances and elastic up to rather high energies. Mislead by this simplicity, one frequently uses 
the leading-order chiral Lagrangian to write the kaon self-energy as = II a LO: (oj) = 

(3a>/4 + Lkn)P/ f 2 , which determines the kaon spectrum as a solution of the Dyson equation, 
ft) 2 — q 2 — m\ — IIs:(ft), q) = 0, with hik being the kaon mass, / ~ 92 MeV standing for the pion 
decay constant and the kaon-nucleon E-term £^,y ~ 350 MeV. However, in Ref. [18] it was shown 
that if the kaon self-energy is evaluated in terms of the real parts of the s- and p-wave kaon-nucleon 
scattering amplitudes obtained from the non-perturbative solution of the Bethe-Salpeter scattering 
equation [19] with the kernel given by the chiral perturbation theory, then the kaon self-energy can 
be parameterized as n^ NP) (ft),^) = (l.lm^ — ft)+ 0.2^ 2 /m^)46.8MeVp/po . We see the scalar 
and vector terms in I I^'° : and I I A XI> have opposite signs. Thus, the chiral Lagrangian does not 
correctly describe the kaon-nucleon scattering process if evaluated in perturbation theory. 

The world of antikaon-nucleon scattering is much richer: it is a coupled-channel system as 
the KN scattering is inelastic already on the threshold due to the coupling to Kf\ and TtZ channels. 
There is a bunch of resonances, among which the s-wave A(1405) resonance dominates the KN 
scattering in the isospin-0 channel close to the threshold. In the p-wave, the resonance £(1385) 
determines the isospin-1 scattering amplitude. The KN interaction is strongly attractive right below 
the threshold this would lead in medium to a shift of the K spectral density for q = 0 towards lower 
energies. There, kaonic modes can couple strongly to the hyperon-nucleon-hole modes (HN 1 ) [2] 
built by ground-state hyperons A(1116) and £(1195). The complicated interplay of anti-kaons and 
hyperon-resonance propagation in medium asks for self-consistent calculations [20, 21]. There are 
attempts to implement the anti-kaon spectral functions in transport calculations [22]. As argued 
in Refs. [1,2] anti-kaons leave the fireball created in a HIC at the last stage of it evolution, when 
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the baryon density is rather low 0.5po < pg < Po- Therefore instead of the spectral function one 
frequently uses an effective scalar potential changing the K mass as m* K = m k + U^Pis /Po with the 
parameter Ug = —(70-150) MeV. 

3. Minimal statistical model for strange particles 

We assume that in HIC a thermalized nuclear system (a fireball) is formed. This hot and 
dense system participate in a hydrodynamical expansion, which lasts until the moment of freeze- 
out characterized by the values p# f Q and 7f 0 . Henceforth, in-medium particle thermal momentum 
distributions become distributions of free-streaming particles. To simplify, we assume the fireball 
to be spatially uniform and characterized by a time-dependent temperature T(t), baryon density 
pri(t) and volume V(t). Because of the high production thresholds the strange particles are most 
efficiently produced at the early hot and dense stage of the fireball evolution. At SIS-SPS energies 
the fireball is baryon rich, and antistrange particles (kaons, K + and K°, mesons) have longer mean 
free paths than strange particles (antikaons and hyperons). Therefore kaons can easily move off 
the production point, and either leave the fireball immediately or first thermalize via elastic kaon- 
nucleon scattering and then leave at some intermediate stage. Because the population of strange 
particles is very small, even if the kaon stays in the fireball for a while, there is little chance that it 
meets antikaons or hyperon and is absorbed by them. Thus, in the course of collision the amount 
of strangeness of the fireball grows. The accumulated strangeness is redistributed among various 
hadronic species, like K~, K°, A, £, Z, £2 and their resonances, which are released at the fireball 
breakup. (At highest RHIC and LHC energies the strange antibaryons will also contribute.) Thus, 
the kaon yield measured experimentally or calculated theoretically in some hadro-chemical models 
can be used to normalize the abundance of strange particles. 

We will demonstrate how this approach works in several examples given below. 

3.1 Strangeness at the SPS energy scan. K + /n + horn 

The excitation function of the ratio (K + )/(n + ) which exhibits a sharp peak (a “horn”) at 
projectile energies of 30 AGeV was intriguing experimental result obtained in the energy scan 
performed at the CERN SPS [23]. Statistical models were not able to reproduce the observed 
shaipness of the peak [24] unless the very broad a meson resonance was included [25]. The 
excitation functions of ( K + )/ (n *"), (K )/ (n~) and (A)/ (n) are difficult to reproduce in hadronic 
transport codes either. In Ref. [26] we calculate relative abundances of strange particles in HIC 
using a hadronic kinetic model for strangeness production. Being interested only in the ratios of 
total yields it is enough to study the evolution of the spatially averaged densities of individual 
species only. The density of kaon is then determined by the equaiton 

d df =PK (~V^c) + ^ gain "^ loss ’ ^ gain = ZlT8rr pi ‘ pj ’ (3 ‘ 1} 

where V is the volume of the system, t is the time in the co-moving frame. The first term on 
the right-hand-side is for the density change due to the fireball expansion. The term 2% gam is the 
kaon production rate with the summation running over non-strange hadron species with densities 
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Figure 2: Excitation functions of the strange particle ratios calculated in [26] in comparison with the 
data [23, 29], 


Pi and the ( vpOjj) standing for the reaction cross section averaged with the relative velocity of 
the colliding particles over the relativistic Boltzmann distributions. The most important reaction 
channels with nB, BB (B = N.A) and MM (M = n.p) entrance channels are included. The last 
term, &i oss , on the right hand side of (3.1) represents the rate of kaon annihilation processes. As 
demonstrated in [26] for AGS-SPS energies, B^i oss <C & ga i n and can be neglected. Equation (3.1) is 
integrated starting from some initial kaon density due to primordial kaon production in collisions of 
incident nucleons. It was estimated from a compilation of kaon production data in nucleon-nucleon 
collisions [27]. 

For the evolution of the energy density and the baryon density we use parameterizations con¬ 
sisting of an initial accelerating period and a later scaling expansion with power-law dependence of 
the density on time, whereby the evolution of the baryon density determines actually the expansion 
of the fireball volume. By exploring a range of parameters we can vary the evolution between the 
well-known Bjorken and Landau hydrodynamic solutions of the fireball expansion. 

The input parameters are the total fireball lifetime and the initial energy density. They, together 
with the final state densities, determine the whole evolution scenario. The energy density and 
number densities in the final state are chosen so that they correspond to the values extracted within 
the statistical model [28] from the chemical freeze-out fit of the experimental data. 

Solving Eq. (3.1) we obtain the total amount of anti-strangeness (S=+l) produced in the fireball 
during its expansion. Species with negative strangeness must balance the total strangeness of the 
system to zero. Reactions which just swap the strange quark between them are quick. Therefore, 
we assume that these species are in chemical equilibrium with respect to each other, while the 
strangeness-weighted sum of their densities is given by the requirement of strangeness neutrality. 

The resulting particle ratios are depicted in Fig. 2. Our hadronic kinetic model [26] is able to 
describe satisfactorily the excitation function of the ratios K + /n + , K~/n~, and A/ti. However, 
the ratio Z /A remain underestimated in the entire range of collision energies. This can indicate, 
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perhaps, that the £ baryons decouple from the fireball at an earlier stage of the fireball evolution at 
an higher temperature. In Ref. [26] we have not included in-medium effects neither for hyperons 
nor for mesons, which are not important at the fireball temperature characteristic for this energy 
range. 

3.2 <j) production 

The study of 0 meson production in HICs provides complementary information on collision 
dynamics and in particular on the evolution of the strange subsystem. Interactions of 0 mesons, 
consisting mainly of s and s quarks, with non-strange hadrons are suppressed due to the Okubo- 
Zweig-Iizuka rule. Hence, amplitudes of OZI-forbidden reactions are typically orders of magnitude 
smaller than those of OZI-allowed ones. Since the OZI suppression weakens the 0 production only 
in the ordinary hadronic matter and would be lifted in the quark-gluon medium, the strong, order of 
magnitude, enhancement of the 0 yield was proposed in [30] as a signal of the quark-gluon plasma 
formation. Experiments at AGS and SPS energies had indeed indicated an enhanced 0 yield albeit 
to a lesser degree [31]. In [34] the enhancement was explained by contributions from the OZI- 
allowed process with strangeness coalescence KK —> 0p and KA —y <j)N, provided the 0 meson 
mass is strongly reduced in medium. Surprisingly strong enhancement of the 0 yield was observed 
at the beam energies about 2 GeV per nucleon [32]. Such a large 0 abundance cannot be described 
by the transport model [33] where 0s are produced in OZI-forbidden reactions BB —y BB(j) and 
kB —t 05 (5 = N, A) with the dominant contribution from pion-nucleon reactions. Note that the 
strangeness coalescence process could not contribute much at these energies, since kaons have a 
long mean free path and most likely leave the fireball right after they are created. 

In [35] we proposed a new mechanism of the 0 production on strange particles, nY —> 07, 
KN —> 07 with 7 = A, £. In contrast to the coalescence reaction, here the strangeness does not 
disappear from the system and plays the role of a catalyser. The presence of K mesons is unneces¬ 
sary. It was found [35] that the catalytic reactions can be operative if the maximal temperature in 
nucleus-nucleus collisions is larger than 130 MeV and the collision time is larger than 10 fm. 

If the 0 production in HIC is indeed promoted by other strange particles one could expect 
correlations among them to be seen in observables. The centrality dependence of the 0 production 
was argued in [35] could be such a one. The number of produced 0 mesons can be estimated 
as N<j, Cl conv Kp -fflcatNpp 3 , where N pp is the mean number of projectile participants and the 
first term corresponds to conventional processes like kN N(j), and the second term is due to 
the catalytic reactions. The extra power N pp appears in this term because the number of strange 
particles in the fireball scales with the fireball expansion time to which in turn is proportional 
to the system size / oc y 1 / 3 « Npp 3 for the case of ideal hydrodynamics. For the experimental 
ratios measured in [36] for Au+Au collisions at I 1.7 AGeV we find N$/N K ~ an pp + bn pp , and 
No /N k i ~ a! + b n pp , where n pp = N pp /A, a,a' and b,b' parameterize the relative strength of 
conventional and catalytic processes, and A is the number of nucleons in the colliding nuclei. Fits 
of the experimental data by these two expressions are presented in Fig. 3 (left and middle panels). 
We see that the best fits are reached when both parameters are activated (dash-dotted curves). 

A link between 0 meson and hyperons could be also seen in particle rapidity distributions. 
If, following Ref. [37], we assume that the rapidity distributions of particles stay unchanged after 
some initial stage when nuclei are passing through each other, then the rapidity distribution of 
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Figure 3: Observables in 0 production indicating the link among 0 mesons and hyperons. Left and middle 
panels: centrality dependence of ratios tj)/n and (!) / K v . Data points are from [36] correspond to Au+Au 
collisions at 11.7AGeV. Right panel: width of the (!) rapidity distributions in PbPb collisions versus the 
beam rapidity [37]. Curves show the rapidity widths calculated from rapidity widths of reactant species. 

0 mesons produced in reaction 1 + 2 —y 0 + X can be approximated as the product of rapidity 
distributions of colliding particle species 1 and 2. In Fig. 3 (right panel) we show the experimental 
root-mean-square of the rapidity distributions of 0s produced in PbPb collisions [37] in comparison 
with the rapidity distributions width for 0s stemming from K'K — y 0, TtA —y 0A and K A — y N(j) 
reactions. We observe that reactions involving hyperons produce the distributions width in better 
agreement with the experimental results. 

4. Rare strangeness: S puzzle at GSI energies 

We turn now to HIC at energies close to the strangeness production threshold, like those stud¬ 
ied by FOPI and HADES collaborations at GSI. Now the strangeness production is really a very 
rare event and we have to respect the strangeness conservation on event by event basis. Depending 
on the number of s quarks remaining in the fireball (recall, s quarks leave the fireball in form of K 
mesons at early stage of the collision), various strange hadrons and their combinations can be ob¬ 
served in the final state. Hence the chemical equilibrium conditions are different for the events with 
different numbers of produced kaons. To handle this complication we divide the totality of events 
with strangeness production into classes of events with one, two, three and so on, ss quark pairs 
created. We call them n-kaon events and denote the probability of creation of exactly n strange 
quark pairs as The probability is determined by the Poisson distribution = W" e w /n\, 
with the integral probability of the creation of one ss pair given by IF = A V f ' where Vf 0 is the 
fireball freeze-out volume determined by the freeze-out baryon density pg f G and the overlap vol¬ 
ume of the initial colliding nuclei. So, for a symmetrical collision of two nuclei with the atomic 
number A we can write Vf 0 (b) = 2A F (b / b milx ) /pax, , where Fix) is the overlap function param¬ 
eterized, e.g., in Ref. [38]. The parameter A can be fixed by the total multiplicity of K mesons 
observed in an inclusive collision, .M K + = (nP s - )/(1 + r \). Here the parameter rj = (A — Z) jZ 
quantifies the braking of the isospin symmetry in HIC. (For example in collisions Ar+KCl studied 
by HADES T] ~ 1.14.) The anglular brackets stand for the averaging over collision impact param¬ 
eter (...) =2 /q max dbb (... )/b^ aax , where the integration over the impact parameter runs up to the 
maximal possible value b max . As a result for the parameter A we find A = (1 + rj ) /(V)^ 3 ). 
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The statistical probability that strangeness will be released at freeze-out in a hadron of type a 
with the mass m a is given by the standard Gibbs’ formula 

» 2 r ft 

Hr, 


Pa=z s s a V foPa =z s s °V fo v a e Ba ^ f(m a ,Tf 0 ), f(m,T ) = 




m 


2n 2 


(4.1) 


where B u is the baryon number of the hadron, the degeneracy factor v a is determined by the 
hadron’s spin I a and isospin G a as v a = (2I a + 1) (2G a + 1), and Ko is the MacDonald func¬ 
tion. The baryon chemical potential at freeze-out is determined by pn.to — — 2f 0 In (4 [/(mjv, Pf 0 ) + 
4/(m A , 7f 0 )] /p B ,fo) , where my is the nucleon mass. A isobars are also treated as stable particles 
with the mass m\ = 1232 MeV, and the small contributions of heavier resonances, hyperons, and 
anti-particles are neglected. 

The quantity zs in Eq. (4.1) is a normalization factor. It is related to the probability of finding 
one s quark in the hadron a. The factor z s s a follows from the requirement that the sum of prob¬ 
abilities for production of different strange species and their combinations, which are allowed in 
the finale state, is equal to one. The factor z s s a depends on how many strange quarks are produced. 
Hence, it is different in single-, double- and triple-kaon events. Therefore, we introduce the nota¬ 
tion Pa"' 1 = (zs ] y a Vfo Pa , where the superscript n indicates to which class of events this probability 
and the zs factor belong: 

(i) In a single-kaon event one 5-quark can be released as K, A or £. Hence, the normaliza¬ 
tion condition for the probabilities (4.1) reads 1 Vf 0 (pp + p,\ + p £ ) = 1. The multiplicity Ma 1 ’ 
of strange hadrons of type a = {K .A.Y) produced in such single-kaon events is given then by 


Vfo p a ■ The isospin factor g a takes into account the asymmetry in the yields of 


particles with various isospin projections induced by the global isospin asymmetry of the collision 
r]. For t] = I this factor reduces to the standard 1 / (2 G a + 1). 

(ii) In a double-kaon event there can be one S baryon and all possible combinations of kaon and 
hyperon pairs. For double-kaon events the normalization condition is determined by equation 
z s 2 V£,(Pk +Pa + Pz) 2 +z$ Vf 0 pz = 1. The factors 2, which would appear here at the cross 
terms after opening the brackets, reflect the number of combinations with which two s quarks can 
be released as a given combination of hadrons, for example, KA, KY and EA. Now the multi- 


( 2 ) A. A. Z, 

plicity of the particle a with one s quark produced given by M a ’ =g a 2py>prL b p;r where 
a,b = {k\ A.£}. We take here into account that the hadron a can be produced in pair or in various 

combinations with other strange hadrons. In both cases the quark combinatoric factor 2 has to be 

( 2 ) ( 2 ) ( 2 ) 

included. The multiplicity of produced S baryons is M~ = gv PY PA> . 

(iii) In the very rare triple-kaon event, when three K ' mesons are produced, the normalization con¬ 
dition is 4 3)3 y fo (PR + Pa + /h) 3 + 3 z [ P 3 VI{pk + Pa + pz)pz + z^VfoPa = F Factors 3 hidden 
in this relation show in how many different ways three s quarks can be distributed between three 
hadrons (e.g., KAA, KKA) or two hadrons (e.g., ZK, EA). The multiplicity of hadrons with one s 
quark is M® = g a 3 P^P^ [( PP +P {3) - 1 - p(3) ^ 2 - 1 - p(3) 
m 1 3) = g E 3P s (3) P| 3) (P^ 3) + P| 3) + P^ 3) 

Having the normalization factors and the chemical potential at our disposal, we can calculate 
the multiplicity ratios of strange particles as functions of the freeze-out density and temperature, 


,( 2 ) d ( 2 ), 


j(2) 


P/C + Py ) +P E ]■ For hadrons with two s quarks we find 
and for the Q baryon = P^p P^\ 
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Figure 4: The strange particle ratios (dashed lines) calculated with Eq. (4.2) in comparison with the available 
experimental values (symbols). 
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These relations are applied in Ref. [39] to the description of multiplicity ratios of strange particles 
detected by the HADES collaboration in Ar+KCl collisions at 1.76AGeV [40]. The main effect 
from the exact strangeness conservation prove to comprise a reduction of the ratio R~,/a/k+ by 
factor ~ 0.5, the ratio Rq/a/k- /k+ by ~ 0.2 and the ratio Rq/s/k+ by factor ~ 0.4 compared to the 
results obtained without exact strangeness conservation in each event. The inclusion of in-medium 
potentials for hyperons and K mesons, as discussed in Sect. 2 is important, leading to much better 
agreement with the experiment. For the K potential Uf> — —15 MeV the freeze-out baryon density 
Pb, to = 0-6po and temperature 7f 0 = 69 MeV are found, which agree well with the parameters 
obtained from the analysis of the pion production in the same energy range [11], The resulting 
particle ratios are depicted in Fig. 4 in comparison with the experimental ratios. We see that the 
calculations are in good agreement with the experiment for the rations of singly-strange particles. 
The multiplicity of S baryons is, however, strongly underestimated. 

Several possible sources of the E enhancement were discussed in Ref. [39] with the conclusion 
that to get any substantial increase in the number of E’s we have to assume that these baryons are 
not absorbed after being produced and their number is determined by the rate of direct production 
reactions, as, e.g., for dileptons. There two types of reactions: the strangeness creation reactions, 
like KN —> KZ, nH —> KZ, and strangeness recombination reactions, like KH —> Zn and HH' —> 
ZN. The former ones have very high thresholds and, therefore, are not operative at SIS energies. 
Strangeness recombination reactions are secondary processes involving two strange particles. They 
can be promoted statistically. We argue that once several s quarks appear in hadronic system its 
more energetically favorable to “store” them together in a multistrange state, E or D. In Fig. 5 we 
show the mass spectrum of strangeness -2 and -3 states in channels with different baryon numbers 
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Figure 5: Mass spectrum (energy E) of strangeness -2 states (left panel) and -3 (right panel) with baryon 
numbers 0, 1, and 2 counted from the non-strange ground states with the energy Eq (shown on top of the 
figure). 


B. We see that in both single and double baryon channels the states with the E and Q baryon are at 
the bottom of the spectrum. This means that Ss and Ds would play a role of a strangeness reservoir, 
being filled with a decrease of the temperature. The only sink is the reaction Z/V —y AA which has, 
however, a relatively small cross section [41]. 

5. Multistrageness at LHC 

We turn now the the highest energies reached in HICs so far. The ALICE collaboration at 
CERN LHC reported [42] the momenta distributions of particles produced in Pb+Pb collisions 
at s/snn = 2.76ATeV. We analyze these data with the help of the DRAGON Monte Carlo event 
generator [43]. The DRAGON based on the blast wave model [44] but includes not only the 
particles stable with respect to strong interactions, but also hadron resonancnes, handling altogether 
277 hadron species. The particle distribution in the transverse momentum (pr) and rapidity y space 
is given by the integral drN/dprdy = / d 4 xS(x,p) over the source function 

S(x,p) d A x = di m> —— n±(pflU? - Pi)/T)B{ 1 - r/R)rdrd(p S(t - T 0 ) Tdtd rj . (5.1) 

Assuming that the fireball expands dominantly in a longitudinal direction we use here the lon¬ 
gitudinal proper time T = \Jt 2 — z 2 and space-time rapidity r/ = jln((f + z)/(t — z)). A position 
in the transverse plain is defined by polar coordinates r, (p. The statistical distribution function 
n±(jt) = I/(exp (A) ± 1) is used with sign + for fermions and — for bosons, whereas di stands 
for the spin degeneracy of the hadron species. All members of isospin multiplets are treated 
separately. The source function (5.1) assumes sharp freeze-out along the hypersurface t = To 
and uniform density distribution within the radius R. This implies that the freeze-out time does 
not depend on the radial coordinate. The fireball expansion is parameterized by the velocity 
field = (cosh r\ t cosh rj, sinhrj, cos ( 0 ,sinhrj r sin <p, cosh rj, sinh 77 ) where the transverse veloc¬ 
ity is such that v t = tanh rj, = rj f(r/R) n . In this relation rjj defines the transverse flow gradient 
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Figure 6: Transverse momentum distribution of n , K ± , p (left panel) and A, E and <4 (right panel) in 
Pb+Pb collisions at s/snn = 2.76ATeV in comparison with the experimental data [42] 


and n selects various transverse velocity profiles. The mean transverse velocity is then equal to 
(v t ) = 2r\f/(n + 2). The transverse size R and the freeze-out proper time To influence total normal¬ 
izations of transverse momentum spectra. However, in this study we ignore those and hence we 
have no sensitivity to these geometric parameters. 

With the help of DRAGON we fitted spectra of pions, charged kaons, protons and anti-protons, 
together with spectra of K° mesons and As for the most central collisions [42]. The relative particle 
abundances in the collision are fixed at the moment of the chemical freeze-out characterized by the 
temperature T c h = 152 MeV and the baryon chemical potential /i« ~ 0 [45]. The successive decay 
of hadronic resonances is explicitly treated in the code. The set of parameters (T, T]f, n ) in Eq. (5.1) 
controlling the stage of the kinematic freeze-out of the fireball is optimised for the best agreement 
with the data. The details of the fitting procedure can be found in [46]. 

For the most central collisions it was found T = 98 MeV, t]j = 0.88, and n = 0.69, that cor¬ 
responds to the average transverse velocity (v t ) = 0.65. The quality of fits is illustrated in Fig. 6. 
On the left panel we show the spectra of most abundant particles p, n ~, K ± , which are, in gen¬ 
eral, very well described by the model in the wide range of transverse momenta. Only the pion 
data demonstrate an enhancement at low p t , which cannot be fully accounted for by contributions 
from decaying resonances. This deficit, we believe, can be resolved by the introduction of non¬ 
equilibrium potential for pions. This would naturally occur if the hadron gas chemically freezes 
out at a temperature around 150-160 MeV and then cools down while keeping the effective ratios 
of individual species constant. We estimated that the pion chemical potential at kinetic freeze-out 
temperature might reach values around 100 MeV. This is not enough for Bose-Einstein condensa¬ 
tion but modifies the spectrum considerably. 

The spectrum of As is also well reproduced as seen in Fig. (6), right panel. However for 
multi strange particles, which were not included in the common fit, the fit results predict spectra 
which are in worse agreement with the data. The experimental data are in general steeper than the 
predictions. The data would thus indicate lower transverse expansion velocity and earlier freeze-out 
of multistrange baryons. 
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6. Conclusions 

In the present talk we discussed several topics related to the production of particles with several 
strange quarks. 

(i) For AGS-SPS energies we show that using a hadronic kinetic model, which describes 
the K ' meson production, we arc able to reproduce excitation functions of the multiplicity ra¬ 
tios ^/n ± and A/n. In the same time, the E/A ratio remains systematically underestimated for 
AGS-SPS energies. 

(ii) The enhanced yields of 0 mesons at SIS and AGS energies, as we argue, could ask for 
account of a new type of reactions in which the 0 production is assisted (catalyzed) by a hyperon 
or a kaon. Such a correlation among 0s and strange particles could manifest itself in the centrality 
dependence of the 0 yield and in the 0 rapidity distribution. 

(iii) In collisions with a small number of produced strange particle (SIS energies) the exact 
strangeness conservation in each collision event must be explicitly preserved. This implies for 
instance that S baryons can be released only in events where two or more kaons arc produced. This 
has an influence on the observed mean number of Es. The observed enhancement of E production 
may mean continuous production of multistrange hyperons when they freeze out right after they 
arc produced in the system. 

(iv) For collisions at the LHC energies we reconstruct a final state of with the help of DRAGON 
Monte Carlo generater, based on the blast wave model supplemented with the decay of hadronic 
resonances. It is found that pion, kaon, (anti)proton and lambda spectra can be satisfactorily de¬ 
scribed by the common freeze-out temperature and flow parameters. However multistrage baryons 
seem to decouple at earlier stage of the collision at higher temperatures and weaker transverse flow. 
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